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Abstract.  We  have  recently  observed  several  features  from  a  randomly 
rough  dielectric  film  on  a  reflecting  metal  substrate  including  a  change  in 
the  spectrum  of  light  at  the  satellite  peaks,  the  high-order  correlation, 
and  enhanced  backscattering  from  the  grazing  angle.  In  this  paper,  we 
focus  on  the  enhanced  backscattering  phenomena.  The  backscattering 
signal  at  small  grazing  angles  is  very  important  for  vehicle  re-entrance 
and  subsurface  radar  sensing  applications.  Recently,  we  performed  an 
experimental  study  of  far-field  scattering  at  small  grazing  angles,  espe¬ 
cially  enhanced  backscattering  at  grazing  angles.  For  a  randomly  weak, 
rough  dielectric  film  on  a  reflecting  metal  substrate,  a  much  larger  en¬ 
hanced  backscattering  peak  is  measured.  Experimental  results  are  com¬ 
pared  with  theoretical  predictions  based  on  a  two-scale  surface  rough- 
ness  scattering  model.  ©  2004  Society  of  Photo-Optical  instrumentation  Engineers. 
[DOI:  10.1117/1.1646410] 

Subject  terms:  enhanced  backscattering  effect;  coherence  effect;  rough  surface 
scattering. 

Paper  030252  received  May  27,  2003;  revised  manuscript  received  Sep.  18, 
2003;  accepted  for  publication  Sep.  18,  2003. 


1  Introduction 

Interference  effects  with  diffuse  light  have  been  studied  for 
a  long  time.  A  recorded  observation  of  the  phenomenon 
was  made  by  Newton  about  three  centuries  ago,1  with  a 
description  of  the  appearance  of  a  series  of  colored  rings 
that  occurs  when  a  beam  of  sunlight  falls  on  a  concave, 
dusty  back-silvered  spherical  mirror.  The  phenomenon  was 
explained  later  by  Young,2  Herschel,3  and  Stokes4  by  con¬ 
sidering  the  interference  between  two  streams  of  light:  one 
scattered  on  entering  the  glass  and  the  other  scattered  on 
emerging  from  the  glass.  The  attention  to  this  phenomena 
was  attracted  once  and  again  many  times  during  the  last 
two  centuries  (see,  for  example,  Raman  and  Datta,5  and  the 
bibliography  in  de  Witte6  and  Francon7,8).  During  the  last 
several  decades  this  problem  again  has  received  attention 
owing  to  numerous  applications  in  modem  optics  and 
microelectronics9"16  and  radio-physical  methods  of  solving 
remote  sensing  problems  for  natural  layered  media.17"25 

For  scattering  of  light  from  a  rough  dielectric  film  on  a 
reflecting  substrate,  there  are  three  main  kinds  of  trajecto¬ 
ries  that  give  rise  to  (a)  Quetelet  fringes,  (b)  Selenyi 
fringes,  and  (c)  enhanced  backscattering.  A  typical  Quetelet 
ring  pattern  consists  of  a  series  of  elongated,  colored  dif¬ 
fuse  rings.  The  white  ring,  corresponding  to  the  zero  order 
of  interference,  passes  through  both  specular  and  back- 
scattering  directions,  and  as  the  angle  of  incidence  is 
changed,  new  colors  emerge  from  the  center. 

Besides  the  Quetelet-type  ring,  there  is  another  kind  of 
interference  effect  that  has  been  observed  in  the  scattering 
of  light  from  dusty,  back-silvered  mirrors.  The  Selenyi 
rings26  present  a  different  kind  of  behavior  under  oblique 
illumination;  there  is  no  zero-order  ring,  and  the  rings  are 
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always  centered  around  the  normal  to  the  sample.  The  oc¬ 
currence  of  such  rings  has  been  explained  in  terms  of  the 
interference  between  waves  scattered  back  directly  from 
the  top  of  the  scattering  layer  without  entering  it,  and 
waves  reflected  by  the  mirror  after  first  having  been  scat¬ 
tered  when  entering  the  film. 

One  of  the  most  interesting  phenomena  associated  with 
the  scattering  of  light  from  a  randomly  rough  surface  is  that 
of  enhanced  backscattering.  This  is  the  presence  of  a  well- 
defined  peak  in  the  retroreflection  direction  in  the  angular 
distribution  of  the  intensity  of  the  incoherent  component  of 
the  light  scattered  from  such  a  surface.  It  results  primarily 
from  the  coherent  interference  of  each  multiply  reflected 
optical  ray  with  its  time-reversed  partner.  For  rougher  me¬ 
tallic  surfaces,  where  the  rms  roughness  is  about  1  to  3  /nn, 
the  enhanced  backscattering  was  measured  with  high- 
sloped  single-scale  Gaussian  random  rough  surfaces  27  For 
weakly  rough  metallic  surfaces,  where  the  rms  roughness  is 
about  0.01  //.m,  the  enhanced  backscattering  was  predicted 
in  a  perturbation  calculation  and  subsequently 
measured.  8,29 

The  scattering  of  light  from  a  one-dimensional  randomly 
rough  dielectric  film  deposited  on  a  flat  reflecting  substrate 
has  been  studied.30,31  In  particular,  here  we  investigate  the 
appearance  of  well-defined  fringes  in  the  angular  distribu¬ 
tion  of  the  diffusely  scattered  intensity  and  their  depen¬ 
dence  on  the  angle  of  incidence,  the  roughness  of  the  film, 
and  the  film’s  mean  thickness.  It  was  found  that,  for  slightly 
rough  films,  the  angle  of  incidence  modulates  the  intensity 
of  the  fringes  but  has  no  effect  on  their  angular  position. 
For  rougher  films  there  is  less  contrast  in  the  pattern,  and 
the  fringes  move  with  the  angle  of  incidence  in  such  a  way 
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that  there  are  always  bright  fringes  in  the  specular  and 
backscattering  directions.  Eventually,  for  very  rough  films, 
the  fringe  pattern  disappears,  and  a  well-defined  back- 
scattering  peak  emerges  in  a  retroreflection  direction. 

The  measurement  of  the  scattering  of  electromagnetic 
waves  from  a  randomly  rough  surface  at  grazing  angles  of 
incidence  presents  a  challenging  problem;32  this  is  due  at 
least  in  part  to  the  fact  that  if,  for  example,  a  one¬ 
dimensional  random  surface  is  illuminated  by  a  beam  of 
finite  width  W  (see  Fig.  1),  its  intercept  with  the  mean 
scattering  surface  A  =  W/cos(0;),  where  is  the  angle  of 
incidence  measured  counterclockwise  from  the  normal  to 
the  mean  scattering  plane,  increases  to  a  very  large  value  as 
6 1  approaches  90  deg.  For  example,  if  0f=  89  deg,  A 
=  57.3W.  We  have  to  select  a  small  beam  size  W — about 
1.5  mm.  Sample  length  L  of  the  random  surface  should  be 
large  enough  to  compromise  the  grazing  angle  edge  effect; 
L  is  therefore  chosen  to  be  200  mm. 

The  theoretical  interpretation  of  observed  interference 
effects  usually  is  based  on  two  main  approaches.  For  large 
roughness,  where  rms  height  exceeds  essentially  the  light 
wavelength,  the  tangent  plane  method,  based  on  Kirchhoff 
diffraction  integral,  is  employed.  For  the  problem  of  scalar 
waves  scattering  by  a  rough  surface  it  was  applied  first  by 
Brekhovskikh33  and  generalized  by  Isakovich34  for  electro¬ 
magnetic  wave  scattering  by  random  rough  surfaces.  The 
detailed  statement  of  this  method  and  its  further  develop¬ 
ment  can  be  found,  e.g.,  in  the  monographs  of  Beckmann 
and  Spizzichino,35  Bass  and  Fuks,36  Ishimaru,37  and 
Voronovich.38 

In  the  opposite  limiting  case  of  smooth  interfaces  with 
small  roughness,  where  rms  height  is  much  smaller  than  the 
light  wavelength,  the  small  perturbation  method  is  em¬ 
ployed.  In  the  papers  of  Elson, 39-42  Elson  et  al., 43-45 
Bousquet  et  al.,46  and  Duparre  et  al 47  it  was  developed  for 
electromagnetic  wave  scattering  (“vector  scattering 
model”)  by  rough  interfaces  of  multilayer  dielectric  stacks. 
This  model  was  successfully  applied  and  developed  by 
Amra48-50  and  Amra  et  al.51-54  for  a  theoretical  interpreta¬ 
tion  of  numerous  experimental  results. 

In  this  paper,  we  report  the  observation  of  enhanced 
backscattering  at  grazing  angles  for  a  very  smooth  dielec¬ 
tric  film  on  a  reflecting  metal  substrate.  For  theoretical  in¬ 
terpretation  of  obtained  experimental  results,  we  apply  the 
combination  of  two  mentioned  above  methods:  the  small 
perturbation  theory  for  surface  microroughness,  and  the 
tangent  plane  method  for  large-scale  roughness,  which  can 


be  considered  as  layer  thickness  variations  as  well.  This 
approach  is  similar  to  the  “two-scale  scattering  model” 
commonly  used  for  interpretation  of  radio  wave  and  sound 
scattering  by  natural  rough  surfaces  (see,  e.g.,  Bass  and 
Fuks36).  The  main  difference  is  that  in  radio-physics  and 
acoustics  applications,  the  large-scale  roughness  effects  in 
the  scattering  process  by  the  angular  modulation  of  the  mi¬ 
croroughness  scattering  indicatrix  are  obtained  in  the 
framework  of  the  small  perturbation  method.  Here,  the 
large-scale  roughness  is  assumed  to  be  very  gentle  with 
negligible  slopes,  and  they  effect  scattering  only  by  phase 
disturbances  of  multiple  reflecting  waves  propagating  in¬ 
side  the  layered  structure.  Section  2  will  introduce  this  the¬ 
oretical  analysis,  followed  by  Section  3,  the  experimental 
results,  and  Section  4,  the  summary. 

2  Theoretical  Analysis 

2.1  Small  Perturbation  Theory 

Here,  we  consider  the  small  perturbation  theory  application 
to  the  specific  case  of  a  plane-stratified  dielectric  medium 
with  only  one  rough  surface  separating  it  from  the  free 
space  (in  particular,  it  can  be  the  “upper  interface”  of  the 
multilayer  dielectric  stack).  For  this  scattering  problem, 
there  is  no  need  to  use  the  general  form  of  vector  scattering 
theory  developed  in  Refs.  39-47  for  multilayer  dielectric 
stacks  with  an  arbitrary  number  of  rough  interfaces.  We  can 
employ  for  this  specific  case  the  simplified  and  more  com¬ 
pact  version  of  this  theory  reported  in  Fuks  and 
Voronovich55  and  Fuks.56 

It  is  assumed  that  the  surface  roughness  height  s(r) 
(where  r={x,y}  is  a  radius- vector  in  the  plane  z  =  0  of  the 
medium  boundary)  is  much  smaller  than  the  incident  wave¬ 
length  X.  The  specific  (related  to  the  unit  area  of  surface 
z  =  0)  scattering  cross-section  cr^ky  ,kz)  of  the  plane  in¬ 
cident  wave  with  wave  vector  k,  (see  Fig.  1)  and  polariza¬ 
tion  state  fi=p,s  polarization  corresponds  to  the  direc¬ 
tion  of  electric  vector,  which  is  perpendicular  to  the  plane 
of  incidence  xOz,  and  p  polarization  corresponds  to  the 
electric  vector  that  lies  in  this  plane]  into  the  scattered 
plane  wave  with  wave  vector  ks 
=  ( k0  sin  0S  cos  <ps  ,Icq  sin  6S  sin  <ps  cos  0S)  and  polarization 
state  a-p,s  [5  polarization  corresponds  to  the  electrical 
vector  that  is  perpendicular  to  the  plane  of  scattering,  which 
forms  the  azimuthal  angle  <p  with  the  plane  of  incidence 
xOz,  and  p  polarization  corresponds  to  the  electrical  vector 
that  lies  in  this  plane]  can  be  written  in  the  following  form 
[see  Eq.  (5.3)  in  Fuks56]: 

o-%{ks,ki)  =  TTkl\s- 1  \2\fafi\2S%(q±),  (1) 

where  k0  =  27rl\  is  the  wave  number,  e  =  e(0)  is  the 
upper-limit  value  of  the  dielectric  permittivity  s(z)  in  the 
medium  (z<0),  q=ks-k(  is  a  vector  of  scattering,  is 
its  projection  on  the  plane  z  =  0,  and  5s(q1)  is  a  spatial 
power  spectrum  of  surface  roughness,  which  can  be  intro¬ 
duced  as  a  Fourier  transformation  of  roughness  auto¬ 
correlation  function  W(p)  =  s(r  +  p)<;(p): 
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//  W(p)eippdp,  (2) 

where  the  bar  .77  denotes  the  statistical  averaging  over  the 
ensemble  of  random  function  s(r),  and  factors  fap  (which 


are  proportional  to  the  amplitude,  or  length  of  scattering) 
are  given  by  the  expressions  (5.6)  to  (5.9)  from  Fuks56: 

fss  =  [  1  +  Rs(  */)][  1  +  Rs(  0S)] cos  <p,  (3) 

fps=- U  -Rp( 0S)J[  1  +Rs( ®i)]o os  0S  sin  <p,  (4) 

fpp  =  - -[  1  +  Rp(  0,)][  1  +  Rp(  ^)]sin  0,-  sin  ^  [1  Rp(  0,)] 
X[l-Rp(6s)]cos  0jcos8scos<p  (5) 

fsp= LI -Rp(0i) J[1  +/?s(0f)]cos  6 1 sin  <p.  (6) 


Here,  Rs  and  Rp  are  the  specular  reflection  coefficients 
from  the  lower  (z<0)  stratified  media  into  the  upper  half¬ 
space  (z>0)  for  s  and  p  polarizations,  correspondingly. 

In  the  particular  case  of  uniform  (homogeneous)  media 
with  e(z)  =  e  =  const,  when  there  is  no  volume  scattering  at 
all  in  the  half-space  z<0,  coefficients  Rp  and  Rs  coincide 
with  the  usual  Fresnel  reflection  coefficients  R0p  and  R0s 
from  the  plane  interface  of  two  homogeneous  media  (e0 
=  1,  z>0)  and  (e,  z>0): 

e  cos  0—  Je  —  sin2  0  cos  0—  Ve  — sin2  0 

°p  e  cos  0+  yfe  —  sin2  0  cos  0+  yfe  —  sin2  0 

(7) 


In  the  general  case  of  an  arbitrary  stratified  medium,  the 
reflection  coefficients  Rp  and  Rs  can  be  represented  in  the 
form 


R 


p,S 


l+Rop.sR'p.s’ 


(8) 


where  Rp  s  are  the  reflection  coefficients  from  the  buried 
layers  (Rp  s  =  0  for  the  uniform  e(z)  =  const  homogeneous 
half-space  z<0). 


•y/e  —  sin2  0—  \fsi~- sin2  0 
11  yfe-s in2  0+  \fep- sin2  0 

(9) 

8j  ye  —  sin2  0—  e  \[e [  —  sin2  0 
lp  Si\[e- sin2  0+S\fe j  — sin2  0 


and  <p(  0)  =  2 k0H\fe  -  sin2  0. 

A  detailed  analysis  of  the  intensity  spatial  distribution 
pattern,  originated  from  scattered-by-roughness  z  =  s(r) 
and  multiply-reflected  (from  planes  z  =  0  and  z=~H) 
wave  interference,  was  done  by  Fuks56:  the  interference 
rings,  angular  positions,  their  polarization  dependence,  pe¬ 
riods  of  intensity  oscillations  as  functions  of  parameters  H 
and  X,  etc.  These  theoretical  results  are  in  good  agreement 
with  experiments  carried  out  with  the  perfect  Fabry-Perot 
parallel-slided  plates  (H= const)  and  the  very  small  surface 
settled  scatterers  (see,  e.g.,  de  Witte6).  But  in  some  experi¬ 
ments  (see  Kaganovskii  et  al.57)  an  essential  disagreement 
with  this  theory  was  discovered.  It  was  shown  that  large- 
scale  roughness  (LSR)  can  be  the  main  reason  of  destruc¬ 
tion  of  the  interference  between  some  types  of  waves, 
which  leads  to  the  “survival”  of  only  one  specific  set  of 
interference  maxima  and  the  suppression  of  others.  The  the¬ 
oretical  analysis  conducted  in  Kaganovskii  et  al.57  was  re¬ 
stricted  to  consideration  of  interference  of  once-reflected 
waves  only.  For  a  very  low  grazing  angle  of  incidence 
rr/2  —  0,  the  multiple  wave  reflections  into  the  resonator 
formed  by  two  planes  z=0  and  z  —  —  H  can  play  the  lead¬ 
ing  role  in  forming  the  scattered  field  intensity  spatial  dis¬ 
tribution,  and  in  particular,  in  forming  the  backscattering 
intensity  peak. 

To  investigate  the  effect  of  LSR  on  the  scattered  inten¬ 
sity  distribution,  we  assume  that  the  successive  wave  re¬ 
flections  inside  the  layer  (0>z> -//)  take  place  every 
time  from  a  horizontal  plane,  as  depicted  in  Fig.  1,  without 
changing  the  reflective  angles  0\  and  0's ,  correspondingly, 
before  and  after  scattering  by  the  rough  patch  of  the  upper 
boundary.  The  layer  thickness  H  is  different  at  the  different 
points  of  reflection,  as  shown  also  in  Fig.  1.  The  solution  of 
this  modeling  problem  can  be  represented  in  the  form  of 
Eq.  (1)  with  amplitudes  fa/3 ,  which  can  be  obtained  from 
those  given  by  Eqs.  (3)-(6),  with  the  following  formal  pro¬ 
cedure.  In  fa/3 ,  represented  by  the  above-mentioned  equa¬ 
tions,  the  factors  [1  ±RStP\  can  be  rewritten  in  the  following 
form,  using  Eq.  (8): 


2.2  Scattering  in  Presence  of  Large-scale 
Roughness 

Here,  we  apply  the  above  equations  to  the  simplest  case  of 
layered  structure:  the  homogeneous  dielectric  layer  of 
thickness  H  and  permittivity  e  lying  on  the  homogeneous 
half-space  (z^H)  with  a  complex  dielectric  permittivity 
constant  8j.  The  reflection  coefficient  Rp<s( 0)  from  this 
structure  for  every  linear  polarization  (s,p)  is  given  by  Eq. 
(8),  where  Rp  ,( 0)  can  be  written  in  the  form  R'p  s(0) 
=Rlp  s(0)exp[i<p(0)],  where  /?lpii(0)  is  the  Fresnel  reflec¬ 
tion  coefficient  from  the  interface  of  two  media  with  dielec¬ 
tric  permittivity  constants  e  and  £i : 


1  ±R= 


R0±l 

Ro 


1+  *°r_L 

1 +R0R' 


(10) 


Here,  for  brevity,  we  omit  the  arguments  (0,  and  0S)  or 
only  their  subscripts  (i  and  s),  and  the  polarization  sub¬ 
scripts  (p  and  s)  in  reflection  coefficients,  to  avoid  confu¬ 
sion.  We  expand  Eq.  (10)  in  a  series  of  R'=Ri  exp[i<p] 
powers,  which  is  equivalent  to  field  representation  as  a  se¬ 
ries  of  multiple  reflections.  We  substitute  instead  tup,  the 
phase  of  the  wave  re-times  reflected  from  the  undulated 
interface  between  layer  and  substrate  rup=:>'Lnk=-i4>k,  where 
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cpk(  0)  =  2k$Hk yfe  -  sin2  ft  and  {Hk}  is  the  set  of  layer 
thicknesses  in  the  specularly  reflecting  points: 


1  ±R= 


R*±l 


R  n 


l±(*o-i)X  (  RoRiV 

n  =  0 


Xexp  1 2  • 

**=1 


(ID 


Now  we  can  consider  //  as  a  random  function  of  two  vari¬ 
ables  (j t,y),  with  a  given  average  value  ( H },  variance 
and  correlation  length  /#. 

Here,  we  assume  that  there  are  no  losses  inside  the  di¬ 
electric  layer,  i.e.,  Im  e=0,  and  present  results  only  for  the 
limiting  case  of  extremely  strong  variations  of  layer  thick¬ 
ness  a2H ,  when  the  corresponding  Rayleigh  parameter  es¬ 
sentially  exceeds  unity,  i.e., 


<o?f> 

^=Vf-  =  C(0,)C(es), 


(14) 


where  function  C(  6)  is  given  by  the  equation 


(W*)2}  =  4£oOtf(e-sin2  0s)  =  (2kl(THcos8's)2>l,  (12) 

where  6[  and  6S  are  related  by  Snell’s  refraction  law: 
sin  ^'=sin  0slti,  n=yfe  is  a  refraction  index,  and  k\ 
=  nk0.  If  this  inequality  holds,  it  is  possible  to  neglect  the 
averaged  value  of  exponents  in  Eq.  (11);  i.e.,  take 
(exp(z<^))=0  in  all  equations  that  appear  from  Eq.  (11). 

The  result  of  statistical  averaging  of  scattered  light  in¬ 
tensity  angular  distribution  over  the  set  of  random  variables 
{ Hk },  which  we  denote  by  the  corner  brackets  (...), 
strongly  depends  on  the  ratio  of  correlation  length  lH  and 
the  distances  A  between  the  sequential  specular  reflecting 
points  (see  Fig.  1).  If  all  distances  A  between  every  two 
arbitrary  reflecting  points  from  the  substrate  exceed  essen¬ 
tially  the  correlation  length  lH ,  then  all  subsequent  specu¬ 
lar  reflections  from  the  lower  layer  boundary  can  be  con¬ 
sidered  as  independent  random  events,  and  consequently, 
the  set  of  <j)k  is  the  set  of  independent  variables.  Thus,  if 
inequality  A >lH  holds,  then  in  all  directions  of  scattering 
given  by  angles  0s,<p,  excluding  the  vicinity  of  back- 
scattering  direction  (0s=0i9  <p-7r),  the  factors  1  ±R(0j) 
and  1  ±R(0S)  are  statistically  independent  and  can  be  av¬ 
eraged  separately.  The  statistical  average  value  of  the  scat¬ 
tering  cross-section  (<r^g)  over  layer  thickness  fluctuations 
SH=H—{H)  is  proportional  to  the  (\fap\2),  which  for  s 
polarization  takes  the  form  [see  Eq.  (3)]: 


C(0)  = 


l  +  r2(l+2r0) 

i-fwi)4 


(15) 


where  r0—  |/?0j|  and  rj  =  |/?1<s| .  In  the  special  case  of  a 
perfectly  conducting  substrate,  when  ^=1,  Eq.  (15)  takes 
the  form 


C(0)  = 


1  -ro(0)’ 


(16) 


and  we  obtain  the  simple  equation  for  the  contrast  coeffi¬ 
cient,  introduced  by  Eq.  (14): 


KSs(6s>6i)  = 


[  1  r0(  #;)][  1  “  ^o(  ^j)] 


(17) 


It  follows  from  this  equation  that  KSS(0S ,  0,-)5=4,  because 
0^r0(^>5)^l,  and  in  the  particular  case  of  e>l,  when 


1  -ro(0)  = 


2  cos  0 


2  cos  6 


Ve-  sin2  ^+cos  0 


we  obtain  for  Kss(  6S ,  : 

e 


Kss(0sA)  = 


cos  0S  cos  0t 


> 1 . 


(18) 


(19) 


(l/«|2)  =  (il+^(^)|2)(|l+^(^)|2)cos2<p.  (13) 

Compare  the  scattering  cross-section  (cr^),  averaged  over 
the  layer  thickness  variations,  with  the  corresponding  value 
ass  for  the  homogeneous  half-space  bounded  by  the  same 
rough  surface.  The  explicit  equation  for  arss  is  given  by  Eq. 
(1)  with  fss  from  Eq.  (3),  where  the  reflection  coefficients 
Rs  have  to  be  substituted  by  R0s ;  i.e.,  if  we  insert  Rfs  =  0 
into  all  the  above  equations.  For  the  ratio  of  these  two 
cross-sections  (which  can  be  named,  according  to 
Kalmykov  et  al.,58  as  a  contrast  coefficient  Kss ),  we  obtain: 


It  is  seen  that  the  average  brightness  of  the  interference 
pattern  due  to  the  substrate  can  essentially  exceed  the  one 
for  a  homogeneous  half-space  (without  substrate)  even  for 
very  strong  variations  of  layer  thickness  SH  when  all  the 
interference  maxima  are  completely  smoothed. 

2.3  Backscattering  Enhancement 

The  backscattering  case  has  to  be  considered  separately  be¬ 
cause  the  dashed  and  dotted  ray  trajectories  in  Fig.  2  are 
fully  congruent  in  this  case,  and  it  is  impossible  to  carry  out 
the  averaging  over  SH  separately  for  each  of  them.  When 
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Os-Oi  and  <p-  it,  and  all  specular  reflecting  points  for 
dashed  and  dotted  rays  coincide,  we  have  to  carry  out  the 
following  averaging: 


ing  angle  or  large  layer  dielectric  permittivity  e,  when  in¬ 
equality  yjs  —  sin2  0,-^cos  6(  holds,  from  Eq.  (25)  it  follows 
that: 


(20) 

where  1  +/^(0Z)  is  represented  as  the  sum  Eq.  (11)  of  in¬ 
dependent  specular  reflections  from  the  undulated  sub¬ 
strate. 

Omitting  the  bulky  derivations,  we  present  here  only  the 
final  expression  for  the  contrast  coefficient  K0ss  in  a  back- 
scattering  direction  as  a  result  of  averaging  in  the  limiting 
case  of  very  strong  layer  thickness  fluctuations  SH ,  when 
the  inequality  Eq.  (12)  holds: 

K0ss= ^  ={(1  +fi)[(l  +  r\){  \+A)  +  2>r\r  0] 

+  2rj[(l  -A)2+2A(3-A)]}(1-A)-3,  (21) 

where  A  =  (r0ri)2  and  0,  is  supposed  to  be  an  argument  for 
all  reflection  coefficients.  Compare  K0ss  given  by  Eq.  (21) 
with  the  indicatrix  contrast  coefficient  Kss  given  by  Eq. 
(14)  for  directions  of  scattering  9S  close  to  backscattering; 
i.e.,  inserting  0,=  0,  and  k,=  -k/5  we  can  estimate  the 
excess  y  of  the  backscattering  peak  of  {a°ss(  -  k, ,  -  k,*)  over 
the  surrounding  background: 


Kqss 


(22) 


Equation  (21)  takes  the  simpler  form  for  the  specific 
case  of  a  perfectly  conducting  substrate  when  rx  =  1 : 


^Oss 


2(3 -r0) 
(l-z-o)3- 


(23) 


Compare  this  expression  with  the  indicatrix  contrast  coef¬ 
ficient  Kss  given  by  Eq.  (17)  for  directions  of  scattering  9S 
close  to  backscattering;  i.e.,  inserting  0S  =  9V : 


4 

"“[l-r0(  */)]*' 


(24) 


We  can  estimate  the  excess  y  of  a  backscattering  peak  of 
(o-®s(— k,-,  —  kj)  over  the  surrounding  background,  in  this 
specific  case: 


(3-rp) 

r_2(l-r0)’ 


(25) 


It  is  easy  to  see  that  the  general  phenomenon  of  back- 
scattering  enhancement,  in  the  specific  problem  under  con¬ 
sideration,  appears  as  a  backscattering  peak  enhancement, 
which  for  (l-r0)<§l  can  essentially  (many  times)  exceed 
the  background,  in  contrast  to  the  well-known  volume  scat¬ 
tering  problem,  where  this  enhancement  can  achieve  only 
the  value  of  several  units.  In  the  specific  case  of  low  graz- 


y= 


yje  —  sin2  9 1 
2  cos  9 1 


>1. 


(26) 


It  is  worth  emphasizing  that  none  of  the  above  equations 
for  contrast  coefficients  Kss  and  K0ss ,  as  well  as  for  back- 
scattering  peak  enhancement  y,  depends  on  the  statistical 
parameters  of  upper  surface  roughness,  and  in  particular,  on 
its  spatial  power  spectrum  5s(qx).  For  a  nonabsorptive  di¬ 
electric  layer  with  lme=0  neither  do  they  depend  on  the 
mean  layer  thickness  ( H )  and  its  variance  a2H ,  if  the  in¬ 
equality  Eq.  (12)  holds. 


3  Experimental  Results 

The  fully  automated  bidirectional  reflectometer  shown  in 
Fig.  3  is  used  to  measure  the  fraction  of  incident  light  re¬ 
flected  by  the  sample  into  incremental  angles  over  its  field 
of  view.  It  uses  illumination  from  laser  sources  at  0.633  fim 
and  enables  measurements  to  be  taken  for  any  combination 
of  incident  and  reflected  angles  over  the  entire  plane,  ex¬ 
cept  for  a  small  angle  (about  0.5  deg  away  from  the  retrore- 
flection  direction)  in  which  the  source  and  detector  mirrors 
interfere.  A  laser  beam  passes  through  a  polarizer  and  is 
interrupted  by  a  chopper  and  a  half-wavelength  plate, 
which  enables  rotation  of  the  polarization  of  the  beam. 
Then  it  is  directed  toward  the  sample  by  a  folded  beam 
system  that  collimates  it  into  a  parallel  beam  up  to  25  mm 
in  diameter.  For  our  measurement,  the  beam  size  is  set  to 
1.5  mm.  The  sample  is  viewed  by  a  movable  off-axis  pa¬ 
raboloid  that  projects  the  light  reflected  by  the  sample  onto 
the  detector  via  a  polarizer  and  a  folding  mirror.  Four  dif¬ 
ferent  polarization  combinations  of  input  and  receiving 
beams  are  recorded.  The  reference  standard  used  for  these 
experiments  is  Lab  Sphere  Gold,  and  the  relative  bidirec¬ 
tional  reflectance  is  measured.  The  signal  is  recorded  and 
digitized  at  each  angular  setting  of  interest  throughout  the 
angular  range  by  an  ITHACO  lock-in  amplifier  and  the  data 
are  stored  in  the  memory  of  a  personal  computer  (PC).  The 
sample  and  the  receiving  telescope  arm  are  separately 
mounted  on  two  rotational  stages  run  by  two  independent 
stepper  motors  that  are  controlled  by  the  PC  via  a  two-axis 
driver. 

The  sample  we  used  is  made  of  smooth  aluminum  that 
was  coated  with  a  dielectric  film  for  high  performance  and 
protection.  The  thickness  of  the  layer  is  approximately  H 
=  5.2  fim.  The  complex  permittivity  ex  of  A1  at  X 
=  0.6328  /xm  is  ex  =  -56.52+21.25/.  The  refractive  index 
of  the  film  is  n=  1.64  (the  dielectric  constant  is  e  =  2.69). 
The  sample  was  characterized  with  a  Dektak  3030  stylus 
machine.  The  stylus  radius  was  approximately  0.15  /mi, 
and  a  total  of  2000  points  over  a  2-mm  profile  were  taken 
for  each  trace.  From  a  statistical  analysis  of  10  profiles  it 
was  determined  that  the  rms  height  of  the  roughness  of  the 
film  was  about  0.06  fi m  and  the  approximate  value  of  Me 
correlation  length  of  the  roughness  was  about  0.3  /mi.  The 
illuminating  source  in  the  experiment  is  a  15-mW  He-Ne 
laser  with  \  =  0.6328  /an.  Since  the  dielectric  film  is 


Optical  Engineering,  Vol.  43  No.  3,  March  2004  563 


Gu,  Fuks,  and  Ciftan:  Grazing  angle  enhanced  backscattering  . . . 


smooth,  most  of  the  energy  goes  in  a  specular  direction, 
thus  a  sensitive  photomultiplier  is  used. 

In  the  experiment,  the  beam  size  is  set  to  W=  1.5  mm. 
Since  the  beam  size  is  small,  the  average  far-field  speckle 
size  is  large.  We  have  to  average  about  100  measurements 
to  obtain  the  far-field  scattering  at  each  angle  by  scanning 
very  small  yaw  and  pitch  directions  of  the  sample. 

Figure  4(a)  shows  the  experimental  results  for  p  polar¬ 
ization,  and  Fig.  4(b)  for  s  polarization  with  Bi 
=  -  89  deg.  There  is  a  large  enhanced  backscattering  peak 
at  near  grazing  angle  05=89deg.  The  ratio  of  the  back- 
scattering  enhancement  peak  over  the  surrounding  back¬ 
ground  at  0y=89  deg  for  p  polarization  is  about  19.6,  and 
the  ratio  of  the  backscattering  enhancement  peak  over  the 
surrounding  background  at  0^=89  deg  for  s  polarization  is 
about  20.4.  However,  in  the  experiment,  the  backscattering 
peak  has  a  finite  width,  and  considerable  uncertainty  ap¬ 
pears  in  determining  the  peak  excess  over  the  surrounding 
background. 

The  dependence  of  K0ss  on  the  angle  of  incidence  0Z-  is 
calculated  for  layer  parameters,  corresponding  to  the  ex¬ 
periment  described  above,  and  shown  in  Fig.  5.  It  is  seen 
that  for  the  very  low  grazing  angle,  the  backscattering  con¬ 
trast  K0ss  can  achieve  values  of  several  thousands.  The  plot 
of  y(0i)  presented  in  Fig.  6  shows  the  layer  parameters 
corresponding  to  the  experimental  curve  in  Fig.  4(b).  For 
0*  =  89  deg,  the  backscattering  peak  excess  y=20  approxi¬ 
mately  coincides  with  the  value  observed  in  the  experi¬ 
ment. 

4  Summary 

We  have  observed  several  features  from  a  randomly  rough 
dielectric  film  on  a  reflecting  metal  substrate  including  a 
change  in  the  spectrum  of  light  at  the  satellite  peaks,  the 
high-order  correlation,  and  enhanced  backscattering  from 
the  grazing  angle.  In  this  paper  we  have  focused  on  the 
enhanced  backscattering  phenomena  at  the  grazing  angle. 


Recently,  the  coherence  theory  predicted  that  the  corre¬ 
lation  in  the  fluctuations  of  a  source  distribution  can  cause 
frequency  changes  in  the  spectrum  of  the  emitted  radiation, 
even  when  the  source  is  at  rest  relative  to  the  observer.59 
Because  of  the  close  analogy  that  exists  between  the  pro¬ 
cesses  of  radiation  and  scattering,  a  similar  effect  is  ex¬ 
pected  to  occur  when  radiation  is  scattered  by  a  static  me¬ 
dium  or  rough  target.  This  phenomenon  also  has  a 
significant  relationship  to  enhanced  backscattering,  which 
is  one  of  the  most  interesting  phenomena  associated  with 
the  scattering  of  light  from  a  randomly  rough  surface.  This 
is  described  as  the  presence  of  a  well-defined  peak  in  the 
retroreflection  direction  in  the  angular  distribution  of  the 
intensity  of  the  incoherent  component  of  the  incident  beam 
scattered  from  such  a  surface.  For  metallic  optical  rough 
surfaces,  the  high-slope  nature  of  the  surface  leads  to  mul¬ 
tiple  scattering  and  this  contribution  enhances  the  scattered 
light.27  For  a  weak  corrugated  metallic  surface,  it  is  due  to 
the  multiple  scattering  of  the  surface  polaritons  excited  by 
the  incident  light  passing  through  the  hills  and  valleys  on 
the  surface  before  they  are  converted  back  into  volume 
waves  in  the  vacuum,  propagating  away  from  the 
surface.28,29  The  coherent  interference  of  such  a  scattering 
sequence  and  its  time-reversed  partner  leads  the  enhanced 
backscattering  in  the  retro-reflection  direction. 

In  our  case,  the  multipath  propagation  is  caused  not  by 
multiple  scattering  by  surface  roughness,  but  the  multiple 
specular  reflections  of  waves  between  two  surfaces:  the  me¬ 
tallic  substrate  and  the  upper  boundary  of  the  dielectric 
film.  To  obtain  the  effect  of  backscattering  enhancement,  it 
is  enough  to  consider  only  single  scattering  by  a  dielectric 
film  micro-roughness  and  perform  statistical  averaging  over 
the  large-scale  roughness,  with  negligible  slopes,  which  is 
equivalent  to  the  film  thickness  smooth  variations. 

The  scattering  of  electromagnetic  waves  from  randomly 
rough  surfaces  near  the  grazing  incident  angle  is  of  wide¬ 
spread  interest  in  various  physical  situations,  in  particular 
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Dielectric  Film  on  Metal  Substrate 


(a) 

Dielectric  Film  on  Metal  Substrate 


(b) 

Fig.  4  Experimental  results  for  (a)  p  polarization  and  (b)  s  polariza¬ 
tion. 

coastal  and  shipborne  radar  backscatter.  It  is,  however,  a 
very  difficult  problem,  with  issues  such  as  sea  spikes  and 
unexpectedly  large  parallel-parallel  scattering  cross- 
sections,  still  awaiting  a  satisfactory  analytical  solution. 


Fig.  5  Backscattering  contrast  coefficient  Koss . 


Fig.  6  Backscattering  peak  enhancement  y. 


Backscattering  signals  at  small  grazing  angles  are  very 
important  for  vehicle  re-entrance  and  lidar  signature  appli¬ 
cations.  For  a  randomly  weak,  rough  dielectric  film  on  a 
reflecting  metal  substrate,  a  much  larger  enhanced  back- 
scattering  at  0S  =  89  deg  is  measured,  which  is  compared 
with  a  theoretical  calculation.  Due  to  Quetelet’s  rings,  the 
energy  of  diffusion  is  redistributed  and  a  large  portion  of 
energy  is  attracted  to  the  retro-reflection  direction  at  the 
grazing  angle.  For  this  reason,  a  large  backscattering  peak 
appears  on  the  grazing  angle. 
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